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Abstract 

The dynamics of the tip of the falhng chain is analyzed. Results of laboratory experiments are 
presented and compared with results of numerical simulations. Time dependences of the velocity 
and the acceleration of the chain tip for a number of different initial conformations of the chain 
are determined. A simple analytical model of the system is also considered. 
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I. INTRODUCTION 



The problem of bodies falling in a gravitational field is so old that hardly anything new 
could be expected to add to it. However, the development of the numerical simulation meth- 
ods has opened for the analysis a few interesting cases, difficult to analyze with analytical 
means. The dynamics of a falling chain is among them. 

There are a few variations of the problem. For instance, one can consider here the chain 
initially gathered in a compact heap located on a table, close to one edge. The motion 
starts when one of the chain ends is brought over the edge. If one assumes that the chain 
leaves the heap without friction, the problem becomes tractable analytically and, as one can 
demonstrate, it falls with a constant acceleration. The surprise is that the acceleration is 
not g, as one would expect, but 

In the variation of the falling chain problem considered here the chain is initially attached 
at both ends to a horizontal support. Then, as one of the ends is released, the chain starts 
falling. The case when the initial horizontal distance Ax between the ends is small in the 
comparison to the chain length L has been considered before^. We repeated the experiments 
demonstrating that, in this case, the end of the chain falls with an increasing acceleration, 
always greater than g\ this apparently paradoxical result can be explained via analytical 
calculations^. In the present paper, we provide the detailed analysis of the problem and we 
extend the study to the case of large Ax. 

What happens when the initial horizontal separation Ax of the chain ends increases, 
in particular when it reaches its largest value Ax = L, was not known. We describe in 
what follows a series of laboratory experiments performed at different horizontal separations 
between the ends of the chain and compare the measurements with the results of numerical 
simulations. The not studied before case, in which the initial distance between the ends of 
the chain equals L (i.e. when the chain is initially maximally stretched) proves to be very 
interesting. 

II. THE FALL OF THE TIGHTLY FOLDED CHAIN ANALYTICAL SOLUTION 

To get an intuitive insight into what we can expect in the experiments with the falling 
chain, let us consider first the case Ax = (figure [T]). We introduce here an analytically 
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tractable simplified model, in which we assume that the conformations explored by the 
falling chain can be always seen as consisting of two pieces: a) the falling piece that shortens 
with time; b) the almost motionless piece that elongates with time. Such a division of the 
chain is possible when the initial horizontal separation of the chain ends equals zero and the 
chain consists of infinitely many, infinitely thin segments. In this limit the chain can be seen 
as a folded perfectly floppy and infinitely thin continuous filament. 

Initially both ends of the chain are attached to a point of support O, the vertical position 
of which y = 0. Then, at time t = 0, one of the ends of the chain is released and the chain 
starts falling. Figure Q presents schematically the geometry of the system. 

The dynamics of this simple model can be solved analytically by applying the law of 
energy conservation^.. We assume that the chain has total length L and that its mass M is 
distributed uniformly along it. To simplify final analysis of the results, we introduce a new 
variable h describing the distance of the freely falling tip to its initial position. The h-axis 
is thus oriented downwards, in the direction of the gravitational field. In what follows, we 
shall refer to h as the fall distance. In terms of h, the length and mass of the falling part 
are given by: 

r /, N L-h (L- h)M 

Laih) = ^^, M,{h)= ^ . (1) 

The length and mass of the motionless part are: 

r /,N L + h (L + h)M 

Lb{h) = ^-, M,ih)= ^ (2) 

and the corresponding vertical positions of their centers of mass are given by: 

La{h) 3h + L Lb{h) h + L 
Vcaih) = -h — = yc,ih) = — = —. (3) 

The potential energy of the falling part of the chain expressed in terms of the fall distance 
h, relative to the point y = —L/2, equals: 

rr^,^ ..n^ f n^ L\ Mg{L^ -ALh + 3h^) 

Uaih) = Ma{h)g [y^Sh) + ^ J = ^ (4) 

Analogously, the potential energy of the motionless part of the chain is given by: 

U,{h) = M,{h)g {y^Sh) + ^) = ^^M^ZZ^) . (5) 
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The total potential energy of the chain is thus given by the formula: 

Note that, according to the chosen reference, U{L) = (i.e. the potential energy is zero 
when the tip of the chain reaches its lowest position). 

The kinetic energy of the falling part of the chain is given by: 

_ Maih)v,ihf _ MjL - h)v,{hf 

Assuming that the free end of the chain is initially located at /i = 0, we may formulate the 
law of energy conservation for the falling chain (remember that = 0): 

f/a(0) + f/6(0) = Ua{h) + Ta{h) + U,{h), (8) 

which in view of (jH), (0), and ((Tj) gives: 

1 ^ M{L-h){gL-gh + v,{hf) 

4 ^ 4L ■ ^ ' 

After straightforward simplifications, we find the formulae describing the velocity Vc and 
time tc of the chain tip versus the fall distance h: 



h h 



' 

Let us compare the motion of the freely falling chain and that of a compact body. For- 
mulae describing the motion of the compact body can be also derived from the energy 
conservation law: 

g{]^L^=g(^-L + h)i+\v,{h)\ (11) 

Solving equation (fTT|) . we find the velocity and time of the freely falling compact body 
in terms of the fall distance h: 

h 

/ds 
-j= = (12) 



Figures El and 01 present a comparison between the dynamics of the folded chain and the 
compact body calculated for L = 1 m and the gravitational acceleration g = 9.81 m/s^. 

Figure |2h. shows that the velocity of the chain tip is, at any fall distance h, larger than 
the velocity of the compact body. As expected (figure |2t), the acceleration of the falling 



body is constant and equals g while that of the falling chain tip increases with time, always 
exceeding g. 

Comparing the motions of the falling chain and of the compact body one may ask about 
the difference between the times at which they reach the same fall distance h. Figure El 
presents plots providing a clear answer to this question. 

The analytical results obtained above for the simplified model of the falling chain may 
be confronted with results of laboratory experiments performed on a real chain. The results 
obtained for various initial configurations are the subject of the next section ITTTl 

III. LABORATORY EXPERIMENTS 

We are aiming at the experimental study of the falling-chain dynamics and, more specifi- 
cally, at its comparison with the dynamics of a freely falling weight. In order to point out the 
differences in their trajectories, we designed an experimental setup that makes it possible 
to record the simultaneous motions of the two objects. 

The chain consists of stainless-steel identical segments which are made from rods and 
spheres attached to each other (figure The total length of a segment / = (4.46 ± 
0.01)10^^ m and the diameter of the spheres = (3.26 ± 0.01)10^^ m. In addition, we 
determine the minimum radius of curvature Rmm = (4.8 ±0.2) 10"'^ m, which the chain can 
present without loading any elastic energy. We use a chain of length L = 1.022 m, which 
corresponds to n = 229 segments for a total mass M = (2.08 ± 0.01)10"^ kg. 

The chain is tightly attached at one end to a firm support O by means of a thin thread. 
See figure IHl At the other edge located at point P = {xo,yo), the chain ends with a rod (we 
open and remove the last sphere) to which we tie a thin nylon wire (fishing line, diameter 
0.1 10~^ m). The free-falling weight, a sinker (a lead weight used for sinking fishing lines) of 
mass Ai = 10~^ kg, is then attached to the other end of the nylon wire (length about 5 cm). 
We then make the nylon wire hang to two nails and a thin metallic wire (nickel, diameter 
10~^ m) as sketched in figure El The whole system is adjusted so as to insure that the sinker 
and the end of the chain are at the same altitude as the other end of the chain in O {y = 0). 
It can be displaced horizontally in order to choose the initial horizontal separation between 
the two ends of the chain. As the mass of the sinker Ai is about half the total mass of the 
chain M, the system is almost always equilibrated (in addition, the equilibrium is helped 
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by the solid friction in the contact regions of the nylon wire with the nails and the metallic 
wire). Thus, the initial conformation that is formed by the chain after damping of all the 
oscillations, is close to the catenary curve^^. 

Injecting a large electric current (about 1 A) in the nickel wire results in cutting suddenly 
the nylon wire at the point where they cross; the weight and the end of the chain then 
simultaneously start falling freely under the action of gravity. We point out that they both 
fall with a small piece of nylon wire attached to them. However, as the force that pushes 
the wire against the nails vanishes, the friction force vanishes as soon as the wire is cut. In 
addition, during the free fall, the pieces of wire have no effect on the dynamics as the mass 
of nylon is negligible in comparison to the mass of the sinker or chain. 

The falling chain and weight are imaged with the help of a standard CCD video camera 
(Panasonic WV-BP550/G) and the images are recorded on a video cassette recorder. The 
chosen shutter speed (1/4000 s) is enough for obtaining clear images of both the chain and 
sinker (figure IH)). The movies are digitized afterwards by means of a Macintosh computer 
equipped with a frame grabber board (Data Translation DT2255). Further analysis with an 
image-processing software (NIH-Image) makes possible to recover 50 images per second from 
the movies which are initially made from 25 interlaced images per second. The interlacing 
allows thus to double the time resolution but results in a loss in the spatial resolution, which 
is typically of about 4 mm per pixel. 

The positions of both the falling chain tip and the weight at the times ti, i = 0, 1, 2, at 
which consecutive frames were recorded, is determined from the digitized images. To make 
the discussion of the results as simple as possible, the experimentally determined positions 
of the falling object will be given as the vertical h and horizontal w deviations of its current 
x{t) and y{t) coordinates from their initial values {xo,yo)'- 

w(t) = Xq — x(t), 

(13) 

Ht) = yo-y{t), 

In what follows we shall refer to the variables as the vertical, h, and horizontal, w, fall 
distances. According to their definitions, in the initial stages of the falling process both 
of the falling distances are positive. In all experiments yo = 0, while Xq was changing in 
four steps from 1 m to 0.25 m. Note that since the motionless end of the chain is attached 
to point (0,0), the initial horizontal separation of the chain ends Ax = xq. In view of 
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this equality, in the following we shall be denoting the initial separation by xq. Results of 
the laboratory experiments will be confronted in section IVII with results of the numerical 
simulations presented in section M 

IV. THE MODEL OF THE CHAIN AND ITS EQUATIONS OF MOTION 

One can define a few discrete models of the chain; below we present one of them. Its 
equations of motion will be formulated for the case, in which one of the chain ends is attached 
to the fixed support while the other one is free. Let us note that similar models have been 
considered before^*^. 

Consider a chain moving in a gravitational field. Several assumptions will be made to 
simplify the model. First of all, we assume that the chain is constrained to move only 
in the vertical plane denoted by {x,y). A chain of mass M and length L consists of n 
thin cylindrical rods (in the following we shall refer to them as segments) with masses 
rrii = m = M/n, i = l..n, and lengths li = I = L/n, i = l..n. All the segments are rigid and 
cannot be deformed. Consecutive segments are connected by joints with friction. Figure [7| 
shows the geometric representation of our model. 

In order to formulate the equations of motion, generalized coordinates, which rigorously 
determine the state of the system, must be specified. Following our predecessor^', we decide 
to describe the system using angular coordinates indicating the inclination of the consecutive 
segments with respect to the 

The position of the first element is determined by the angle ipi. Similarly, the position 
of the second element is described by the angle (p2- The global conformation of the chain 
in the plane is uniquely expressed by all angles <^i, i = l-.n. The angles are below referred 
to as generalized coordinates of the system. A generalized coordinate (pi indicates an angle 
between the i-th element of a chain and the horizontal axis x. 

The Cartesian coordinates of the i-th mass center {xi,yi) can be written as follows: 



Using the generalized coordinates we shall derive the Lagrange equations of motion. To 




(14) 




i-l 
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start with, we shall consider the energy of the system. The motion of the chain is considered 
as a combination of translational and rotational motions of its segments. Each segment has 
the moment of inertia Jj = l/12m/^, i = l..n, calculated around the axis perpendicular 
to the {x, y) plane and passing through the center of mass of the segment. Taking into 
consideration the relations given in equation (jl4|) . the kinetic energy of the chain is given 
by: 

1 " 

^ = 2 E ('^ + + ^^^') ' (15) 

i=l 

where the dot represents the derivative with respect to the time t. The potential energy 
of the 2-th segment is given by mgyi, where g is the gravitational acceleration. Thus, the 
potential energy of the chain may be expressed as: 

n 

U = ^mgyi. (16) 

To make our model more general, we introduce damping as a Rayleigh's dissipation 
function^: 

1 " 

^r'^i'^i-^i.if , (17) 

i=l 

where r is the dissipation coefficient. We assume that the joint which connects the first 
element of the chain to the support is free of dissipation. This is equivalent to the assumption 
that (fo = 'fi- Similar definition of dissipation was used by other authora^ii. 

The motion of the falling chain is governed by the system of Lagrange equations of second 
kind: 

+ = 0, 1 = l..n, (18) 



dt \ difi J difi dip 

where C = T — U is the Lagrangian of the system. Applying (fT3j) . (HD), (fTTj) and (fTHj) we 
find the set of n equations describing the motion of a chain: 

n n ^ 

rUijCijipj = - ^ niijSijip] + — (v^i-i - 20^ + ^j+i) - yOiQ, i = l..n, (19) 
i=i j=i ^ 

where q = cos(v9j), = cos{ipi — ipj), Sij = sm(ipi — (pj), ai = n — i + ^ and 

{n-i + l, i = j 

■ 
n - max(z,j) + ^, j 

The next section |V| is dedicated to the results of numerical solving of equation |T^. 



V. NUMERICAL EXPERIMENTS 



Equations of motion derived in the previous section can be integrated numerically 
thus allowing one to simulate numerically the motion of the falling chain. In pres- 
ence of dissipation, the resulting system of equations becomes stiff and requires spe- 
cific numerical methods. We selected the RADAU5 algorithm by Hairer & Wanner 
( http : / /www . unige . ch/~hairer/ software . html ) designed for stiff problems. It is based 
on the implicit Runge-Kutta scheme of order five with the error estimator of order four—. 

We performed a series of numerical simulations aiming to reproduce results of the exper- 
iments described in section IIIII Thus, as the initial configuration of the chain we used the 
discrete catenary curve shown in figure with four different separations between the ends 
of the chain: a) Xo = 1.019 m, b) xq = 0.765 m, c) xq = 0.510 m and d) xo = 0.255 m iden- 
tical with the separations used in the laboratory experiments. Numerical simulations were 
performed with n = 229, L = 1.02 m, M = 0.0208 kg, g = 9.81 m/s^ and time t e [0, 0.5] s. 
The only free parameter left was thus the dissipation parameter r. Varying it we aimed 
at obtaining the best agreement of the numerical results with the laboratory experiments. 
To compare the results, we monitored the distance between the positions of the chain tip 
found in the consecutive frames of the video recordings and the positions found in the nu- 
merical simulations at the same times. The distance between laboratory and numerical data 
obtained in a single experiment is defined as follows: 



are here the horizontal and vertical deviation from the initial position of the chain tip found 
in consecutive frames of the laboratory and numerical experiments, respectively. In order 
to find the optimal value of r providing the best fit for all four experiments a), b), c) and d) 
we determined the total distance 



A depends on the assumed value of r; we have found its values for r in the range from 
10^^ to 10^^. Then, A(r) was analyzed with the use of the least-square algorithm based 
on the procedure SVDFIT-2.. The optimal value of dissipation parameter was found to be 




(20) 



where the denotes the number of analyzed frames. Points {wi, hi) and (iDj, hi) for i = 1..N 



(21) 
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equal r = 2.163 ■ 10 ^ Nms (i.e. for this one, the A(r) reaches its minimum value equal to 
0.02510 m). 

It seems interesting to check how this single value of r fits the data obtained in each of 
the four experiments; table U shows the results. 

In all cases, 6 are relatively small until the chain reach its minimal vertical position (it 
is less then 0.004 m). It becomes much bigger after the chin tip start to raise and it has 
the great influence in the value of A. Figure |H1 presents the experimental vertical h and 
horizontal w fall distances together with their numerical counterparts determined with a 
much smaller time step. 

Consecutive conformations of the falling chain found in the numerical simulations are 
presented in figure El The conformations correspond to the same times at which they were 
recorded in the laboratory experiments. Positions of the falling compact body are also shown 
in the figure. As one can see comparing figures IHl and the shapes of the experimental and 
numerical conformations are almost identical. 

VI. QUANTITATIVE ANALYSIS 

Quantitative analysis of the digital images recorded in the laboratory experiments pro- 
vided us with sets of discrete data representing the vertical, h, and horizontal, w, fall dis- 
tances of the chain tip versus time. As described above, using the data we have found the 
values of the dissipation parameter at which numerical simulations fit best the experimental 
data (tabled). As seen in figure |H1 the agreement is very good. Thus, to analyze the details 
of the falling chain dynamics we shall be using the data obtained at small time steps from 
the numerical simulations. 

First of all, let us analyze the most interesting question of the relation between the time 
dependences of the vertical fall distances of the chain tip and the compact body (figure |H1). 

It is worth noticing that in the case a), where the initial conformation of the chain is 
straight and horizontal, the vertical fall of the chain tip and the fall of the compact body 
are identical up to the moment of time at which having reached its maximum vertical fall 
distance the tip starts moving upwards. That it should be like that becomes clear when 
one notices that during the fall the end part of the chain remains horizontal - its vertical 
motion must be thus identical with the fall of the compact body. Why the end part of 
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the chain remains horizontal is also clear. This happens because the chain displays no 
elasticity and no energy is stored in its bent regions. This phenomenon, noticed in the 
laboratory experiments and confirmed in the numerical simulations, suggests the existence 
of an approximate analytical treatment. So far we have not been able to find it. 

In cases b), c) and d) the vertical fall distance of the chain tip, up to the moment of 
time th^^^ at which vertical fall distance of the chain tip reaches its maximum value hmax, is 
seen to be always ahead the vertical fall distance of the compact body. This observation is 
sometimes summarized by the general statement, that the chain falls faster than a compact 
body. 

The next question that we shall analyze is the time dependencies of the velocity Vc and the 
acceleration of the chain tip. In order to do so we perform a series of numerical experiments 
with xq from a range [0.1, . . . , 1.02) m. All parameters of the numerical simulation are the 
same as defined in previous section. The plot starts at xq = 0.1 m since because of the finite 
length of the chain segments at smaller initial separations the simulated dynamics of the 
chain fall becomes very complex. Similar effects are observed in the laboratory experiments. 

By velocity and acceleration we mean here the moduli of the velocity and acceleration 
vectors. Figure ^1 presents both the variables versus time. As one can see, plots of the 
velocity versus time display distinct peaks. It seems interesting to check how high the peaks 
are (i.e. which the maximum velocities of the chain tips are for different initial conformations 
of the chain) and at which times they are reached. It seems also interesting to ask for which 
initial separation of the chain ends the velocity peak is highest. Answers to these questions 
can be found out analyzing figures ^2 and El 

As well visible in figure ITTh . in accordance with expectations, the peak velocity value 
Vmax becomes largest when the initial separation of the chain ends is smallest (i.e. when the 
chain is maximally folded). On the other hand, contrary to expectations, the velocity peak 
is not smallest at the maximum xq = L initial separation but earlier, at Xq ~ 0.9040 m. 

Figure reveals an interesting fact. The moment of time t^^^^ at which the velocity 
of the chain tip reaches its maximum value precedes in general the moment of time 
at which the chain tip reaches its maximum vertical fall distance h^ax- A reverse rule is 
observed only at the largest initial separations Xq of the chain ends. 

It seems interesting to check how the time th at which the chain tip reaches its lowest 

O "'max i 

position depends on the initial separation xq of the chain ends. The dependence is plotted 
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in figure IT^ revealing that the lowest position of the chain tip is reached fastest when xq ~ 
0.5500 m (i.e. when the initial horizontal distance between the chain ends is approximately 
half of its total length), th^^^ proves to be longest in the case when the chain is initially 
straight. 

The last question we asked analyzing the velocity data was the correlation between the 
value of the peak velocity Vmax and the time at which it is reached. This correlation is 
presented it figure IT^. As we have demonstrated before, the peak of the velocity is highest 
at the smallest initial separation of the chain ends, but one should not draw the conclusion 
that it is thus reached in the shortest time. As seen in the figure the initial separation of the 
chain ends at which the velocity peak is reached fastest amounts to about Xq ~ 0.7000 m. 
See also figure IT^. 

Now, let us analyze the behavior of the acceleration f figure (TUb). As in the case of the 
velocity plots, we also observe here clear peaks. This time, however, they are distinctly 
different in their height. Figure ITTb demonstrates clearly, that the highest peak in acceler- 
ation is also observed at the smallest initial separation of the chain end. Its value, at the 
experimentally studied case of xq = 0.255 m amounts to 7352 m/s^, thus it is about 40 times 
larger than the value observed at Xq = 0.765 m (186.3 m/s^). That such large values of the 
acceleration are realistic was demonstrated by Krehl et. al."^ who studied the dynamics of 
the cracking whip. Figure ITTb demonstrates that accelerations at the lowest positions of the 
chain tip are not the maximum ones. 

It seems interesting to check the relation between the time t^^^^ at which the chain 
tip reaches its maximum velocity and the time at which it reaches its maximum fall 
distance. Figure [T^ presents the relation. As one can see, except for the largest values of 
the initial separation xq the maximum velocity is reached before the chain tip reaches its 
maximum fall distance. 

Figure 112b presents an answer to a differently posed question. We ask about the relation 
between the maximum fall distance of the chain tip hmax and the fall distance h^^^^ at 
which the tip reaches its maximum velocity. As before one can clearly see that in general 
the maximum velocity is reached before the chain reaches its maximum fall distance. 

Parametric relations between a) and Vmax, b) hmax and Vmax found at given values 
of the initial separation xq are plotted in figures IT^ and 113b . As well visible in the figures 
the range of large Xq proves to be very interesting. Small changes of Xq lead here to large 
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changes of t^^^^ and hmax- This range of xq needs a special attention. 

VII. SUMMARY AND DISCUSSION 

A chain falhng in the gravitational field can be seen as a model of other systems such as 
the cracking whip ^°i^^ . At the first sight the above statement may seem not true, since in 
the cracking whip problem, the gravitational forces are in general neglected. Let us notice, 
however, that the end of the folded whip attached to the whip handle is subject to a strong 
acceleration. Changing the laboratory reference frame to the non-inertial frame moving with 
the end of the handle, we introduce into the system strong inertial forces equivalent to the 
gravitational ones. This explains the validity of the initial remark. 

Experiments we performed revealed a few new, interesting facts concerning the dynamics 
of the falling chain. Let us summarize them. 

1. Both the velocity and acceleration dependencies on time display distinct peaks the 
height of which depends on the initial separation of the chain ends. The highest peaks 
are observed for smallest initial separation. There exists an approximate, analytical 
description of tightly folded chain dynamics explaining the origin of the rapid increase 
of the velocity and acceleration. The theory is however unable to predict the finite 
height of the peaks. (In the analytical model both the velocity and the acceleration 
diverge.) 

2. The velocity peak is observed to be reached fastest for initial separation xq = 0.6863 L, 
where L is the length of the chain, whereas its amplitude is smallest for xq = 0.8863 L. 

3. It seems very interesting that in the case in which the initial separation of the chain 
ends is largest, the dynamics of the vertical fall of the chain tip proves to be identical 
with the dynamics of the fall of a compact body. That is should be the case becomes 
obvious when one notices that the end part of the chain remains horizontal during the 
fall. This observation suggests the existence of an approximate analytical treatment. 
It is not known yet. 

4. As a rule, the time at which the chain tip reaches its maximum velocity generally 
comes before the time at which it reaches its lowest vertical position. Only at the 
initial separation of the chain ends larger than 0.8863 L to 0.9608 L the rule is reversed. 
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5. The ratio between the the largest and smallest acceleration peaks is about 166.5, which 
is unexpectedly large. This may have some practical implications since at the times, 
when the acceleration reaches its highest values, forces acting on the chain tip also 
become very large what may lead to a damage of the chain. 

Dynamics of the falling chain hides certainly a few more interesting details. The same, 
even to a larger extent, concerns the dynamics of the falling rope. In the latter case the dissi- 
pation plays a much more important role and elasticity becomes a crucial factor. Laboratory 
and numerical experiments are waiting to be carried out. 
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formed in the experiments with chains built from a finite number of segments is not a catenary 
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initial conformation consists of two exactly vertical pieces. This was the conformation used by 
Schagerl in his laboratory experimental 
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Tables 



TABLE I: Distance between experimental and numerical results 



Experiment: 


5 [m] 


a) 


0.007672 


b) 


0.006964 


c) 


0.005912 


d) 


0.004552 
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Figure Captions 
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FIG. 1: Geometry of the conformation of the tightly folded chain at time i > 0. The position of 
the freely falling chain is described in terms of h. Part a) of the chain is falling down while part h) 
is motionless; we denote by Ca and Cb their centers of mass. 
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FIG. 2: The velocity of the falhng chain tip Vc and the compact body Vf, versus: a) the fall distance 
h and b) time t. 
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FIG. 3: a) The fall distance he of the tip of the chain and the compact body, hj, versus time t. 
b) The difference between times at which the compact body and the chain tip reach the same fall 
distance h. 
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FIG. 5: Sketch of the experimental setup. 
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FIG. 6: Successive conformations of the falling chain vs. time. The left end of the chain remains 
attached to the frame, while the right end is free to fall due to gravity. In b), c) and d), the white 
straight-lines connect the free falling end to the free falling weight for the last five images before 
the maximum extension of the chain (length L = 1.022 m, time spacing between the successive 
images 1/50 s, initial separation between the chain ends: a) xq = 1.019 m, b) xq = 0.765 m, c) 
xq = 0.510 m, and d) xq = 0.255 m). 



22 




FIG. 7: Model of the chain. 
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FIG. 8: The comparison of the vertical h and horizontal w fall distances of the falling chain tip 
found experimentally (circles) and numerically (solid lines). The parabola of the compact body 
fall is also shown (dotted lines). The initial separation between the chain ends: a) xo = 1.0195 m, 
b) xq = 0.765 m, c) xq = 0.51 m, and d) xq = 0.255 m. 
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FIG. 9: Successive conformations of the falling chain vs. time found in numerical simulations. 
Simulations were performed with n = 229, L = 1.02 m, M = 0.0208 kg, g = 9.81 m/s^ and the 
values of r given in tabled The initial conformations of the chain were discrete catenary curves 
with a) xq = 1.0195 m, b) xq = 0.765 m, c) xq = 0.51 m and d) xq = 0.255 m. Positions of the 
freely falling compact body are shown at the right parts of the figures; dotted lines connect them 
with the respective positions of the tip of the falling chain. 
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FIG. 10: The moduli of the velocity and the acceleration of the falling chain tip - numerical 
approximation of the experimental data. Initial separation of the chain ends: a) xq = 1.0195 m, 
b) Xq = 0.765 m, c) xq = 0.51 m, d) xq = 0.255 m. 
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FIG. 11: Moduli of the velocity a) and acceleration b) of the chain tip versus the initial horizontal 
separation of the chain ends. Vmax and amax are, respectively, the maximum velocity and acceler- 
ation reached by the chain tip during its fall, f and ah^^^ are the velocity and acceleration of 
the chain tip observed at the moment of time at which the tip reaches its lowest position. Picture 
b is plotted in logarithmic scale. Gravitational acceleration g is marked with a dashed line. 
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FIG. 12: a) t^^^^ - the time at which the chain tip reaches its lowest position; ty^^^ - the time at 
which it reaches its maximum velocity Vmax- 

The dashed line represents the time tc{L) 0.386722 
at which the velocity diverges in the analytical model considered in section 2. b) hmax - the largest 
vertical fall distance reached by the chain tip; h^^^^ - the vertical fall distance of the chain tip at 
which it reaches its maximum velocity. 
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FIG. 13: a) Parametric plot of the time ty^^^ at which the velocity of the chain reaches its maximum 
velocity Vmax versus the value of the latter, b) Parametric plot of the maximum vertical fall distance 
hmax versus the maximum velocity Vmax reached by the end tip. Bigger circles indicate the points 
at which the values of the plotting parameter xq (i.e. the initial separation of the chain ends) are 
identical to those applied in laboratory experiments. 
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